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ABSTRACT 


The combinations of nuclear, propulsion, and airframe components, related by 
12 independent variables, that yielded maximum payload are presented for various flight 
cruise points. A unit shield which minimized maintenance and operational problems, 
was designed to limit the flight crew dose to 2. 5 rem per 1000 hours of flight. Crash 
safety and afterheat removal problems were not assessed. Payload fraction increased at 
a declining rate as gross weight increased. The value was 5 percent for a 1/2 -million- 

pound (22. 7xl0 4 -kg) airplane and increased to 19 and 22 percent, respectively, for 1- 
1 4 4 

and l^- million-pound (45. 4x10 - and 56. 7x10 -kg) airplanes designed for Mach 0. 8 at 
36 000 feet (10 970 m). Comparison with the C-5A class of airplanes at equal gross 
weights (7. 3x10 lb; 3. 3x10 kg) indicates superior payloads for the nuclear plane at all 
ranges beyond 5900 miles (~llxl0® m) for one-way trips, or 2700 miles (~5xl0® m) for 
round trips without refueling. The important characteristics of major components are 
described, as well as the trade-off networks occurring in determining the best payloads. 
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PRELIMINARY STUDY OF A SUBSONIC NUCLEAR CRUISE AIRPLANE 


HAVING A HELIUM-COOLED THERMAL REACTOR* 
by John L. Allen, Laurence H. Fishbach, and William C. Strack 

Lewis Research Center 

SUMMARY 

A preliminary design-point study of a subsonic aircraft powered during cruise by a 
helium- cooled thermal nuclear reactor has been conducted. The crew radiation dose of 
2. 5 rem for 1000 hours of flight duty per year is only one-half the yearly dose recom- 
mended by the Federal Radiation Council. A unit shield was used to minimize mainte- 
nance and operational problems usually associated with this type of aircraft proposal. 

Both metal-water and hydride-water shields were considered, as well as various 
packing factors representing the ratio of fabricated to ideal shield weight. Crash safety 
and afterheat removal weight penalties have not been assessed. The maximum helium 
temperature was limited to 2360° R (1311 K), and the nominal maximum reactor wall 
temperature was fixed at 2560° R (1422 K). Twelve independent variables were used to 
relate nuclear, propulsion, and airframe components so that a maximized figure of 
merit would include trade-off effects. 

Results indicate attractive payloads for most speed and altitude design conditions, if 
the gross weight is large enough. Payload fraction increased with gross weight at a de- 
clining rate. The value for a 1/2-millio. -pound (22. 7x10^ kg) airplane was 5 percent and 

1 4 

increased to 19 and 22 percent, respectively, for 1- and 1^- million-pound (45. 4x10 - 
and 56. 7xl0^-kg) airplanes designed for Mach 0. 8 at 36 000 feet (10 970 m). Comparison 
with the C-5A class of airplanes at equal gross weights (~7. 3x10 lb; 3. 3x10 kg) indi- 
cates superior payloads for the nuclear plane for ranges beyond 5900 nautical miles 
(~llxl0 m) for one-way trips or 2700 nautical miles (~5xl0° m) for round trips without 
refueling. The payload for the nuclear aircraft could possibly be increased as much as 
50 percent with shield technology advances which are unique to the nuclear concept. 

Other technology benefits could be mutually used by the two types of airplanes. 

Payloads varied only 7 percent for turbine inlet temperatures between 1800° and 

2100° R (1000 and 1166 K) and li percent for helium pressures between 1600 and 2000 psi 
6 fi 2 ^ 

(11x10 and 13. 8x10 N/m ) because other component designs were allowed to reoptimize 
for maximum payload. Because of the higher efficiency of the turbofan cycle, its pay- 
load was 108 to 16 percent greater than for the turbojet depending on the design flight 
Mach number. 


♦Title unclassified. 


INTRODUCTION 


Nuclear-powered aircraft today, as in the past, promise essentially unlimited en- 
durance or range. The extensive work of the past ANP (aircraft nuclear propulsion) 
program (ref. 1) did not lead to the production of a usable airplane. Other studies such 
as CAMAL (continuously airborne missile launcher and low level weapon system), the 
logistic work of references 2 and 3, continued to examine the mission application. The 
projects underwent a transition from weapon system to component Advanced Development 
Objectives to eventual cancellation. Nevertheless, as suggested in reference 4, periodic 
reevaluation of performance potential is justified because of advancing technology, new 
ideas, and new mission needs. 

Some of the changes in technology and application suggest that nuclear aircraft ap- 
pear more feasible than indicated by earlier efforts (also see ref. 5): 

(1) The ability to build very large aircraft, such as the C-5A, favors the nuclear sys- 
tem. Shield weight, whichis the major weight fraction of the propulsive system, increases 
at less than a linear rate with power and therefore favors the larger gross weights. 

(2) Present mission emphasis is for subsonic endurance in excess of 50 hours and 
for large cargo capacities or space for outsize items. Since supersonic capability is not 
required, the installed reactor power, and hence engine size and shield weight, is much 
less because of the better lift-drag ratio and lower velocity. Also, the temperature 
margin between the reactor fluid and the turbine air is greater since the ram tempera- 
ture rise of supersonic flight is avoided. 

(3) Similarly, the ability to build a variety of efficient fanjet engines (high bypass 
and compressor pressure ratios, if needed) helps reduce installed power and size. 

(4) Knowledge concerning shield design techniques and materials has improved. 
Therefore, a general reduction of radiation hazard may be realized by specifying a unit 
shield and a low dose rate. The ground handling and maintenance would be greatly sim- 
plified. Furthermore, the smaller core associated with the choice of a closed-loop sys- 
tem using a high- capacity heat-transfer system such as helium or liquid metals (ref. 2) 
means that the unit shield can be lighter than for the open-direct-air method. 

The objective herein is to illustrate or identify as a first step the relative impor- 
tance and characteristics of major components and their interrelations and sensitivity 
factors. Some estimate of the expected payload weight fraction, delivery rate, and en- 
gine and aircraft design will also be given. However, because of the complexity of a nu- 
clear propulsion system, definitive evaluation requires detailed knowledge about many 
components which is not presently available. The results, then, will form the basis for 
more refined calculations as better information on components becomes available. The 
scope of this design-point preliminary study has been limited to a unit-shielded helium- 
cooled thermal reactor using turbojet or turbofan propulsion in a subsonic logistic type 


aircraft. The dose rate has been set at 2. 5 millirem per hour at the flight deck, which 
for a 1000-hour tour of flight duty is well below the Federal Radiation Council limit of 
5 rem per year per individual. Aircraft size was varied from 1/2 to l4 million pounds 
(22. 7x10 to 56. 7x10 kg). A computer program was utilized to trade off simultaneously 
the performance-weight advantages of as many as 12 major design parameters in order 
to maximize a selected figure of merit, such as payload. 

The engines were assumed to have a tandem arrangement of the heat exchanger and 
chemical (fossil) fuel burner located between the compressor and turbine and within the 
turbine frontal area. The burner was used during the noncruise portions of the flight. 
Simultaneous chemical-nuclear operation was not evaluated. Safety or afterheat removal 
problems were not considered and any related weight penalties would have to be deducted 
from the payloads given herein. 
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area 

aspect ratio 

one-dimensional area ratio 
ratio of wetted surface to wing area 
total engine stream tube area 
maximum fuselage width 
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fan jet bypass ratio = flow 

gas generator 

wing span 

structural wing span 

constant 

drag coefficient 

coefficient of drag due to lift 

coefficient of profile drag 

coefficient of total drag 

coefficient of friction 

thickness ratio correction factor 
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lift coefficient 

specific heat at constant pressure 
roughness coefficient, 0. 0005 
wing root chord 

dose distance between reactor and flight deck 

diameter of total airflow stream tube 

diameter of gas generator airflow stream tube 

maximum fuselage depth 

change in total aircraft energy during climb 

aircraft specific energy 

thrust 

thrust per pound of total airflow 

installed thrust 

sea- level static thrust 

correlation factor 

altitude 

fuel heating value 
independent variable 
kinetic energy 

proportionality constant between lift and drag due to lift 

wing quarter chord intersect location, percent of l ^ 
lift 

cruise design lift to drag ratio 
length of cargo compartment 
fuselage length 

characteristic length of a component 

flight Mach number 

number of engines 

ultimate load factor, 1. 5x2. 5 
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pressure 


AP 

AP 

AP 
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HX 


PE 
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P He 
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Q n 
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RC 

R„ 
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sh 


°ht 

S N 

®vt 

S w 

T 

AT 

T 

ex 


pressure drop 

helium pressure drop in duct 

helium pressure drop in heat exchanger 

potential energy 

pressure ratio across pump 

helium pressure out of pump 

reactor inlet helium pressure, psia; N/m 

fan pressure ratio 

overall pressure ratio 

heat exchanger air-side pressure ratio 

reactor power, MW 

pump power, MW 

correction for variation of profile drag with lift 

sweepback correction to induced drag 

correction for nonelliptical induced drag 
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dynamic pressure, lb/ft ; N/m 

range, miles; m (also gas constant for air) 

rate of climb, ft/min; m/min 

outer radius of shield, R c = R s ^ when shield is external to cargo bay 

when reactor is within cargo bay, R g ^ = 0 

fuselage wetted surface area 

planform area of horizontal area 

nacelle and pylon wetted area 

planform area of vertical tail 
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wing planform area, ft ; m 
temperature, °R; K 
temperature rise across reactor 
reactor exit helium temperature 
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r HX 
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L R, wall 
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V app 

V 0 
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W/A 
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crew 
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eng 


W I 

W 

W 
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F, TC 
F+tanks 
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f,a 

f ,1 


W f . 

fix 


W. 

W 


fuse 
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w; 

W 


W 

w 


gear 

HX 


HXL 


W T 


W lines 


maximum helium temperature, 2360° R; 1311 K 
heat exchanger exit helium temperature 
pump exit (reactor inlet) temperature 

nominal maximum reactor wall temperature, 2560° R; 1422 K 
maximum heat exchanger wall temperature 

heat exchanger exit air temperature and turbine inlet temperature without 
chemical burning 

time 

maximum wing root thickness 
landing approach velocity, ft/sec; m/sec 
flight velocity, ft/sec; m/sec 
weight (mass), lb; kg 

helium flowrate per unit free flow in reactor 
weight of crew 
weight of engines 
weight of chemical fuel 

weight of chemical fuel used from takeoff to cruise 
weight of chemical fuel and tanks 
available fuselage load 

total load causing bending moment at wing root 

weight of fixed equipment and systems such as instruments, communica- 
tions, controls, hydraulic, electric, APIPs, etc. 

weight of fuselage 

gross weight 

weight at beginning of nuclear cruise 

landing gear weight 

weight of heat exchangers 

weight of heat exchangers plus ducts 

weight of payload 

weight of helium lines (ducts) 
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nacelle and pylon weight 
weight of powerplant 
weight of provisions 
weight of pump 
weight of shield 

weight of shield plus penetration penalty 

weight of shield plus penetration penalty plus reactor 

weight of structure (includes equipment and systems) 

weight of W p ♦ W eng + Whx + W Unes 

tail weight 

additive wing weight to account for chemical fuel containment 

wing weight 

available wing load 

flow rate, lb/sec; kg/sec 

longitudinal distance between the quarter chord intercept at the center line 
and at the mean aerodynamic chord 

ratio of stream static pressure to sea level static pressure 

exponent 

ratio of specific heats 

o 

ratio of stream total to sea level static pressure, Pq/ 2116 lb/ft ; 

Pq/ 10. 13X10 4 N/m 2 

overall engine efficiency 

pump polytropic efficiency 

ratio of stream total to sea level static temperature, Tq/518. 6° R; 

Tq/288. 3 K 

wing sweepback angle 

wing sweepback angle, quarter chord 

wing planform taper ratio 

ratio of stream static to sea level static temperature 

o 

stress, psi; N/m (also flight path angle) 
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T 


thickness to chord ratio 


r w j_ 2/4 win S thickness to chord ratio perpendicular to wing quarter chord line 

cp packing factor, shield weight multiple 

Subscripts: 

a air 

c cruise 

d descent 

des design 

He helium 

HX heat exchanger 

o reference point values 

rea reactor 

s shield 

T total 

Superscript: 

average value between end of takeoff and cruise 

METHOD OF ANALYSIS 
General Requirements and Goal 

This section presents some background concerning the problems and choices made 
for the nuclear system, propulsion system, and mission. 

General arrangement of a nuclear propulsive system. - Figure 1 defines the partic- 
ular closed-loop helium nuclear propulsion system considered. Fundamentally, heat is 
transferred to the working fluid from the surfaces of the reactor fuel elements; the fluid 
is pumped through pipes or ducts to a heat exchanger where air from the compressor of 
the propulsive device is heated prior to driving a turbine. A biological radiation shield, 
shown here as a unit type, protects selected regions. Some special arrangement for 
penetrating the shield with the working fluid ducts and preventing radiation from M seeing 
through” the duct is needed and is shown here as a half-turn helix. 

Although the typical arrangement of a nuclear propulsive system described previ- 
ously seems very simple, the design of the complete system is complicated by the inter- 
faces between major components. In a general sense, a performance-weight map exists 
for a family of designs for each major component. However, before overall vehicle per- 
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Figure 1. - Schematic of closed-loop nuclear propulsion system. 

formance can be estimated, the component design points must be established after allow- 
ing interaction between components. Thus, the selected figure of merit for the complete 
system can be maximized rather than trying to assemble independently optimized com- 
ponents. The component performance-weight maps used in this study are based on un- 
published preliminary Lewis data. 

A logistics airplane was selected for this study. This choice of a logistic rather 
than an endurance type subsonic airplane for performance estimation was arbitrary. 
However, the logistic missions do not require consideration of such a broad range of de- 
sign points, payloads, fuselage sizes, etc. as the endurance type. The logistics air- 
craft may well fulfill some endurance applications, but, in general, it is inherently com- 
promised by the large and heavy cargo fuselage. 

An ideal logistics airplane that has almost unlimited range should also be large 
enough to contain outsize items of military equipment as well as the ability to carry an 
appreciable load. As previously stated in the INTRODUCTION section, the variation of 
shield weight with nuclear power favors large airplanes. - For these reasons, the base- 
point gross weight was selected as 1 million pounds (45. 4x10^ kg) or somewhat larger 
than the C-5A. 
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The choice of a closed-loop helium coolant system for transferring heat from the 
reactor to the air in the propulsive system is compatible with the goal of very low radia- 
tion levels, since helium does not become activated. Furthermore, helium is chemi- 
cally inert, has excellent heat-transfer characteristics at high pressures, but has mi- 
croscopic leakage problems, as well as the explosion hazard of pressurized gas systems. 

After trial calculations, the heat exchanger was fitted within the turbine frontal area 
and between the compressor and turbine rather than duct air to heat exchangers in a 
more remote location. A burner for chemical fuel was included in tandem downstream 
of the exchanger. The resulting arrangement displays shaft and bearing problems only 
slightly greater than those encountered when scroll ducts are fitted to the engine to 
transfer air to and from a heat exchanger in a remote location. 

Since the engine weight penalty is nearly the same, ducting a good heat-transfer 
fluid to the heat exchanger seems more logical because a greater variety of aircraft ar- 
rangements are possible due to the smaller duct sizes. 

Since nuclear operation is contemplated only for the cruise portion, the chemical 
burner was incorporated for emergency operations and for the takeoff, climb to cruise, 
letdown, and landing segments of the flight. The noncruise parts of the flight path have 
been given only cursory study since this report is primarily a design-point study. How- 
ever, the procedures adopted are believed to be better than assigning typical values for 
items such as fossil fuel weight. 

Thus, a subsonic logistic turbofan or turbojet aircraft of 1 million pounds 
(45. 4xl0 4 kg) gross weight powered by a unit-shielded helium-cooled thermal nuclear 
reactor with heat exchangers integral with the engine has been selected as a conceptual 
configuration. The goal is to use such a configuration as a model for reflecting relative 
importance, characteristics, interrelations, and sensitivity factors of major components, 
as well as defining areas of interest for airplane design-point performance and size. 

System figure of merit and synthesis . - For the logistics mission, the figure of 
merit for the virtually unlimited range nuclear aircraft system is payload or, in some 
cases, delivery rate, which is the product of payload and velocity. The system weight 
synthesis is composed of three primary parts: the payload, aircraft structure, and the 
powerplant. The composition of these primary parts is as follows: 

Powerplant: 

Engines, with helium pump 
Heat exchangers 
Helium ducting 
Reactor with pressure vessel 
Shield and penetration penalty 
Chemical fuel supply 

Allowance for moderator, shield, and after- heat cooling system 



Structure: 

Wing 

Fuselage 

Nacelles and pylons 
Tail 

Landing gear 

Fixed equipment and systems 
Payload: 

Cargo 

Crew 

Provisions 


Powerplant 

Turbo machinery . - Station designations are shown on the schematic in figure 1 
which can represent either the turbofan or turbojet engine. For compressor pressure 
ratios less than 8, a single spool turbojet is considered for engine weight estimation. 
Otherwise, two- spool units are employed for either the turbofan (one fan spool, one com- 
pressor spool) or the turbojet. The relations used for estimating the engine weight and 
the values used for component efficiency are detailed in appendix A. 

The tandem arrangement of the coolant-to-air heat exchanger and the chemical fuel 
burner permits operation at a higher turbine inlet temperature than would be possible 
using only nuclear power. Augmented turbine temperature operation is employed for 
takeoff and climb to cruise, and it is discussed subsequently in the Airplane Configura- 
tion and Operating Assumptions section. The power required to pump the helium is re- 
moved as shaft power from each engine. 

Reactor . - The model selected for the nuclear power source was a helium- cooled 
thermal reactor with water moderation using refractory metal fuel elements at a nomi- 
nal maximum metal wall temperature of 2560° R (1422 K). A sample performance map 
of the reactor family is illustrated in figure 2 for the standard conditions used through- 
out the report, unless specified, of a reactor inlet pressure of 1750 psia (12. 06x10 
N/m^) and an outlet temperature of 2360° R (1311 K). Values of the helium weight flow 
per unit flow area at various power levels are plotted against reactor radius (which is 
outside the reflectors) and a correlation parameter f used in the shield weight calcula- 
tion. 

Reactor performance calculations usually require long periods of computer time 
and are not suited for the optimization procedure used herein which may require several 
hundred reactor iterations for a single design point. The curve fits and Taylor expan- 





Correlation factor, f 


Figure 2. - Reactor radius against shield weight correlation factor for various power 
levels and weight flows per unit area. Helium inlet pressure, 1750 pounds per 
square inch 02.06x10** N/m^l; outlet temperature, 2360“ R (1311 K) ; nominal 
maximum reactor wall temperature, 2560° R (1422 K). 

Sion given in appendix B provide a rapid means of estimating reactor performance. For 
this reason, the range of design cases presented is limited to values between 100 and 
about 400 megawatts. 

Shield . - Shield weight is affected not only by the shield geometry and materials, 
but also by the size (and shape) of the reactor core for a given power level and how close 
the core or reflector can be to the first layer of shield material. Only unit shields are 
considered so that the airplane will not be activated to any extent that would limit opera- 
tions or maintenance. 

The unit shield weight calculations were made for a range of reactor power levels 
and f parameters for two general classes of materials. The calculations were for a 
dose rate of 2. 5 millirem per hour at the nominal crew station (which varied with fuse- 
lage length) and included air scattering effects for sea level density. The base point 
shield is characterized as ’’currently buildable” and would be a relatively straightfor- 
ward design using, essentially, various layers of metal and regions of water. The ’’ad- 
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100 200 300 400 500 


Reactor power, MW 

Figure 3. - Shield weight as function of reactor power. 

Dose rate, 2.5 millirems per hour; distance, 130 feet 
(39.6 ml. 

vanced" shield materials would utilize layers of various metal hydrides and water; how- 
ever, the chemical properties of the hydrides require special measures for protecting 
the shield units from water, air, and stress. The shields will be referred to as metal- 
water or hydride-water. The calculations were made using procedures called shield 
synthesis and Monte Carlo (SANE for neutron dose and SAGE for gamma ray dose, refs. 

5 and 6). Examples of the variation of the unpenetrated unit shield weight with power 
level and factor f are shown in figure 3 and can be represented by the equation 

w s = (Cj + C 2 f)^ (1) 

The parameters and C 2 depend on the shield materials, whereas the exponent j3 
varied only slightly from a value of 0. 406 to 0. 419 for the ranges of conditions and ma- 
terials considered (see appendix B). Over the power range of 100 to 500 megawatts, the 
effect of varying the core size parameter f from -1 to -2 was to increase the shield 
weight about 7 percent for the metal-water shield and about 9 percent for the hydride- 
water shield. 

Shield packing factor and duct penetration: Additional considerations involve the 
penetration of the unit shield with the ducts containing the working fluid, helium, and the 
necessity of incorporating some equipment within the shield such as valves, motors, etc. 
For the latter item, the increased spacing between the core or reflectors and pressure 
vessel, or the placing of equipment within the shield material, increases the volume to 
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be covered by shield material. This is referred to herein as a M packing factor," and 
since it is really determined by detailed layout design, a simple percentage increase in 
shield weight is used in this design-point study. 

Selection of a concentric helium duct for this analysis, such as that used in refer- 
ence 3, allows a rather simple shield penetration. As shown schematically in figure 1, 
a half-turn helix within a certain bore diameter protects against direct neutron and 
gamma streaming, if some excess shield material is added beyond the original unit 
shield dimensions. The pressure drop and shield weight penalties associated with the 
dimensions of the half-turn helical duct penetration were optimized in a separate study. 
For the 1- million-pound (453 600-kg) airplane, the outside diameter of the concentric 
penetration duct was approximately 18 inches (0. 46 m), with a bore or spiral diameter 
of about 4 feet (1.22 m). 

Variation of shield weight with dose distance: Another shield weight variation is 
that due to dose distance as fuselage length to width ratio or gross weight is changed and 
dose rate maintained at a constant level. The relation used is 


W 


s, D d 330 - 100 log 1Q D d 


W 


s, 130 


118.6 


( 2 ) 


for D d in feet or 


278.4 - 100 log in D h 

L°- 3 (3) 

118. 6 

where D d equals (130/264) in feet or (39. 6/80. 4) in meters. 

Heat exchanger and ducting . - The conceptual design of the helium to air heat ex- 
changer, which was configured to fit between the compressor and chemical burner, was 
a four-pass cross- counter flow type. The tube bundle array (ref. 7) was constructed of 

O 

four layers of 1/4-inch- (6. 35x10 -m-) outside-diameter tubes with a nominal distance 

between headers of four feet (1. 22 m). The tube-wall thickness, and hence internal tube 
diameter, was determined by a "thick tube" creep- rupture stress calculation based on 
the applied pressure and a tube-wall surface temperature determined by heat transfer. 
The stress- rupture equation for the material N155 for a creep- rupture life of 1000 hours 
at the surface temperature is 
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or 


a = 7000 exp [-0. 00447 (T^ max - 2060° R)] psi ' 

• > 

o = 4825 exp [-0. 008046 (T waU> max - 1144 K)] N/m 2 ^ 


( 4 ) 


where T wa ^ max is the maximum heat exchanger wall temperature. 

The weight of the necessary tube bundle, headers, and inner and outer pressure 
shells was calculated and multiplied by a design stress safety factor of 10/7. The ac- 
companying air and helium pressure drops were also determined. Since the primary 
component of the heat exchanger is the tube bundle, the main factors that determine the 
heat exchanger weight are the stress- rupture relation for the selected lifetime, the 
maximum tube-wall temperature, and the helium pressure. For a tube life expectancy 
of 10 000 hours, the front constant of equation (4) becomes 3800 for a in psi or 2620 for 
a in newtons per square meter. 

The best procedure for manifolding the concentric helium lines from the separate 
engines prior to entering the shield has not been established. Therefore, a schedule of 
concentric duct length per engine was estimated from preliminary layouts as the routing 
distance from the shield to the middle engine. The values are given in table I. 

A procedure for calculating the weight of the ducts was followed that resulted in the 
minimum duct weight for a given pressure drop and gas conditions. The total duct pres- 
sure drop includes terms for a nominal number of valves or bends. 

Helium pump . - The power required to drive the helium pumps was extracted as 
shaft power from the engine cycle as shown in figure 1 and in equation (5): 


_ Pump power 
& Total installed power 


T. 

in 

AT 


(PR) p 


1 


( 5 ) 


TABLE I. - SCHEDULE OF 
DUCT LENGTH 


Gross weight, W g 

ft/engine 

m/ engine 

lb 

kg 

0. 5xl0 6 

0.227X10 6 

30 

9. 15 

. 85 

.386 

42. 5 

12.95 

1.0 

.454 

50 

15.24 

1.25 

. 567 

— 

60 

18.28 
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The pressure ratio PR is for the complete system and composed of the heat ex- 
changer, line, reactor, and shield pressure losses. The efficiency of the pump rjp was 
taken as 0. 8. The weight of the pump and associated drive was estimated by means of 


*\ 

W pump ' ^He log e< PR > + °- 0708 C p, He iT *He 


or (in SI units) 


( 6 ) 


W pump “ °- 1086 »He VV »°M PR > + 01784 C p, He iTi He 


Chemical (fossil) fuel allowance . - The chemical fuel used in flying the flight path to 
and from the cruise point and for emergency range was approximated by the procedure 
given in appendix C. For chemical operation during takeoff and climb, the turbine tem- 
perature was 2160° R (1200 K). 

3 

An emergency off-set range of 500 miles (804x10 m) was specified for all designs. 
This distance was the sum of the power-off glide range plus the remainder at cruise 
power conditions. 

The actual mode of emergency flight may well differ from cruise conditions; how- 
ever, this assumption at least provides a common first-order allowance varying with the 
airplane design points considered. No chemical fuel allowance for helium circulation 
during nonnuclear operations has been made due to the uncertainty at this time of the 
amount of circulation required. 


Airplane Configuration and Operating Assumptions 

The procedures adopted for estimating aircraft performance and weight represent a 
selection of existing methods configured to respond to the primary variables in this 
study; namely, design Mach number, altitude, and gross weight. The main items re- 
quired from these programs are the drag and structural weight, so that the powerplant 
performance and its weight can be determined prior to obtaining the payload. 

The general requirements for the airplane configuration include cargo volume, 
safety, air-drop ability, and landing speed. One version of a configuration that meets 
these needs is depicted in figure 4. 

Some of the basic configuration assumptions were the following: 

(1) The reactor-shield assembly would be mounted in the fuselage near the longitu- 
dinal center of gravity as high as possible (to give an unobstructed cargo bay). 

(2) The landing gear would be completely fuselage mounted. 
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0 10 20 30 40 50 60 70 80 90 

Meters 

Figure 4. - Schematic of nuclear cruise, logistic airplane. Gross weight, 1x10^ pounds; (453 600 kg) 
cruise Mach number, 0.8; cruise altitude, 36 089 feet (10 990 m). 

(3) The longitudinal center of gravity was fixed at 55 percent of the fuselage length 
to give a large dose distance. 

(4) A dorsal fin was used to help offset this rearward center of gravity. 

(5) The wing was mounted above the cargo compartment and with the quarter chord 
of the mean aerodynamic chord alined with the center of gravity. 

(6) The "T" tail was used to increase the tail moment arm and improve cargo han- 
dling. 

(7) The engine pods, generally six, were mounted below the wing so that an exter- 
nally blown, high-lift system could be used if necessary. 

Airplane polar curve . - In order to determine the aircraft polar curve, knowledge 
of the profile drag and drag due to lift is necessary. More extensive details are given in 
appendix D, while sufficient background to accompany the results are given herein. 
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Profile drag: Because of the logistics role, a schedule of fuselage width with gross 
weight was established together with a general fuselage geometry (appendix D, eq. (Dl)). 
The surface area of the fuselage group was compiled by summing the contributions from 
the basic body, the reactor-shield fairing, the sponsons for the landing gear, and the 
dorsal fin (appendix D, eq. (D2)). The fuselage group surface area only deviated from 
its dependence on gross weight due to small changes in reactor- shield size. The re- 
maining surface areas, wing, tail, nacelles, and pylons were determined within the op- 
timization program to give the maximum payload for the given design Mach number, al- 
titude, and gross weight. The maximum wing loading was restricted to a value compat- 
ible with a landing approach at 200 feet per second (61 m/sec) at a lift coefficient of 2. 5. 
The wing thickness to chord ratio was constant at 0. 18 up to Mq = 0. 65 and then de- 
creased for higher Mach numbers (see appendix D, eq. (D3)). The wing sweepback angle 
at the quarter chord was a function of both flight Mach number and lift coefficient (appen- 
dix D, eq. (D4)). An additional wing geometry restriction was a planform taper ratio of 
0. 3. The horizontal and vertical tail areas were 25 and 10 percent, respectively, of the 
wing area. The surface area of the nacelles and pylons was formulated as a function of 
bypass ratio (appendix D, eq. (D5)). 

The friction drag was calculated as the sum of the individual components based on 
the Reynolds number for a characteristic length of each and a roughness allowance of 
0. 0005 per unit wetted area. The method of reference 8 for turbulent compressible 
boundary layer was used. Additional corrections were made to the wing and tail profile 
drag as a function of thickness ratio (ref. 9). 

Drag due to lift: The contribution of the wing to the drag due to lift was a function of 
aspect ratio, taper ratio, and sweepback angle. An arbitrary constant was assigned to 
correct for the variation of profile drag with lift (appendix D, eqs. (D8) and (D9)). 

Structural weight estimation . - Although the methods employed for estimating struc- 
tural weight are based on old construction techniques, such as skin-stringer, the magni- 
tude of the answers agrees reasonably well with current predictions for very large air- 
planes. The structural procedures (appendix E) are functions of the same primary pa- 
rameters used to configure the airplane and to compute the aerodynamics; therefore, the 
relative complexity is about equal for both the aerodynamics and weight procedures. 

The usual limit load factor of 2. 5 was used. 

Basically, the procedure calculates the weights of the wing and fuselage structures 
according to the applied loads after independent routines have been adopted for weight 
items such as fixed equipment and systems, landing gear, nacelles, and tail. The dis- 
tribution of certain weight items is decided in advance. The nacelles and pylons, in- 
cluding engines, heat exchangers, pumps, and helium ducts, are assumed to be located 
within 30 percent of the semispan; whereas the chemical fuel is distributed over 90 per- 
cent of the semispan. The cargo, reactor and total shield (includes penetration altera- 


tion), fixed equipment and systems, and landing gear, are located in the fuselage. The 
detailed equations are given in appendix E. 

Engine size . - The installed thrust level at the design point was calculated by means 
of equation (7), which includes an arbitrary increment for a rate of climb potential of 
300 feet per minute (91. 5 m/min) as an allowance for maneuverability and hot days: 


RC W g 


( 7 ) 


Generally, six engines were specified, unless the resulting sizes were unreasonable. 
Engine size was not adjusted for any uniform takeoff distance. Defining the operational 
limitations of nuclear aircraft is beyond the purpose of this report, and hence, the im- 
portance of runway length cannot be evaluated until such questions as special bases and 
overflying certain land areas have been resolved. However, the takeoff distance varied 
from a high of about 8000 feet (2438 m) for the design point Mq of 0. 4 and 20 000 feet 
(6096 m) in altitude to 5300 feet (1616 m) for an Mq of 0. 8 and 36 089 feet (11 000 m) 
in altitude. 


Optimization Procedure 

The purpose of the optimization routine was to define the combination of component 
performances that gave the maximum figure of merit, usually payload. The complete 
airplane design, including the nuclear system, the thermodynamic cycle analysis of the 
engine, and the aerodynamics and weight of the airplane, was determined by 12 inde- 
pendent variables. All other items were dependent or assigned parameters. The fol- 
lowing independent variables were iterated until the maximum figure of merit was found 
for a selected design point (Mach number, altitude, and gross weight): 

Engine: 

Fan pressure ratio 

Overall compressor pressure ratio 

Bypass ratio 

Heat exchanger air pressure drop 
Turbine inlet temperature 
Nuclear system: 

Helium temperature out of reactor 
Helium temperature out of heat exchanger 
Helium heat exchanger pressure drop 


Helium duct (line) pressure drop 
Reactor helium weight flow rate per unit area 
Airplane: 

Cruise lift coefficient 
Wing aspect ratio 

Thus, the maximum payload obtained was not the result of independently optimizing 
the system components but, rather, the combination of component performance param- 
eters that gave the best payload. For example, given the variation of the weight of some 
component with pressure drop, the lowest component weight would not be chosen because 
the accompanying high pressure drop would require more pump power, more installed 
reactor power, heavier shield, more drag, larger engines, etc. Some other self- 
consistent arrangement of independent and dependent variables could have been chosen. 
Further details and the flow chart for the optimizing procedure are given in appendix F. 


RESULTS AND DISCUSSION 
Explanation of Trade-offs for the Optimized System 

The calculation procedure is arranged so that for a selected gross weight the pri- 
mary items that are free to vary are wing area (and its associated geometry) and the 
magnitude of the installed power. As previously mentioned (see Profile drag, p. 18) the 
fuselage surface area is relatively fixed. Thus, as wing area, and hence its drag, 
changes, the installed power changes to provide the necessary thrust. Because the fig- 
ure of merit is payload, the trade-off is primarily between the weight of the installed 
powerplant (which includes the items listed in System figure of merit and synthesis , 
p. 10) and the wing weight. Changes in fuselage weight are small compared with those 
for the wing; however, for a given size, the fuselage weight does reflect the loads ap- 
plied by the body contents and the wing. For a selected design point and gross weight, 
the wing weight varies primarily with wing loading and aspect ratio. The powerplant 
weight is primarily determined by trading off the shield weight against the weights of the 
engines, heat exchangers, and chemical fuel. The smaller shield weights due to smaller 
reactor cores as produced by higher helium temperature rise and flow rates per unit 
area must be balanced against the effects of both the increased heat exchanger siz,e and 
pumping power extracted from the engine cycle. Iterative corrections are made for a 
number of effects, such as nacelle-pylon drag and weight. The terminology "optimized 
system" as used hereinafter refers to the largest calculated payload that results for the 
combination of independent variables (listed under Optimization Procedure, p. 19) after 


each has been separately perturbed in sequence numerous times (see appendix F). Ex- 
ceptions to this terminology, such as fixing the value of an independent variable, are 
noted in the subsequent results and discussion. The maximum design temperature of the 
reactor wall was maintained at 2560° R (1422 K), the maximum helium temperature was 
2360° R (1311 K), and, unless noted, the helium pressure (pump discharge) was 1750 
psia (12.06X10 6 N/m 2 ). 

Performance of the optimized system . - The performance calculated for an opti- 
mized system using metal-water shields with a packing factor of 1. 33 for a range of 
gross weights from 1/2 to 1-^ million pounds (0.227x10® to 0. 567x10® kg) and a flight 
spectrum from Mach 0. 4 to Mach 0. 8 at altitudes from 20 000 to 45 000 feet (6096 to 
13 720 m) is presented in figures 5 to 7 and tables II and HI. Results for the 1-million- 

4 

pound (45. 4x10 -kg) airplane are shown in figure 5. 

A payload fraction of 22. 8 percent of the gross weight was obtained for a design al- 
titude of 36 089 feet (11 000 m) at Mach 0. 645 with only a small reduction (10 000 lb; 

4536 kg) out to Mach 0. 73 (fig. 5(a)). 

At a design altitude of 20 000 feet (6096 m) the peak payload fraction increased to 
26 percent at Mach 0. 50, and decreased to about 15. 5 percent at 45 000 feet (13 720 m) 
and Mach 0. 675. These variations in payload are due to changes in powerplant and 
structural weight fractions (see tables n and HI). These numerical results should be 
viewed as "gross" payload, inasmuch as various undefined penalties, such as crash 
hazard protection devices, are likely. 

Converting these payload fractions to delivery rate (defined as the product of cargo 
weight and velocity), expressed as millions of ton- miles per hour (as in fig. 5(b)), shows 
that peak delivery rates occur at higher Mach numbers than peak payload, particularly 
for the 20 000-foot (6096-m) design case. The comparative peak magnitudes show an in- 
crease of 10 percent over the 36 089-foot (11 000- m) value for the 20 000-foot (6096 m) 
altitude, and a decrease of 32 percent for the 45 000-foot (13 720-m) design. 

The wing loadings (at cruise) after the takeoff and climb chemical fuel has been con- 
sumed, the lift-drag ratios, and total engine stream tube area associated with the pre- 
vious spectrum of design cases are shown in figure 5(c). (Wing loading is just another 
way of representing the independent variable, cruise lift coefficient. ) 

The results represent the balance between powerplant and structural weight frac- 
tions chosen by the optimizing procedure for maximum payload, as explained in the pre- 
ceding section (Explanation of Trade-offs for the Optimizing System). The wing size 
was designed by the landing condition only for the 20 000-foot (6096-m) case at Mach 
numbers greater than about 0. 525. For the other two altitude cases, the effects of 
sweepback angle on wing weight was offset by reverting to a lower wing loading and as- 
pect ratio (see tables n and in) for Mach numbers greater than about 0. 7. 

This trend of decreased payload and optimum wing loadings between Mach 0. 7 and 


TABLE n. - MAJOR WEIGHT FRACTIONS AND GENERAL DESIGN PARAMETERS OF OPTIMIZED SYSTEM 
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Design flight Mach number, Mq 


(c) Wing loading, engine area, and lift-drag ratio. 
Figure 5. - Concluded. 


0. 8 was investigated at the design altitude of 36 089 feet (11 000 m). The payload de- 
crease (0. 22 to 0. 19) was mainly the result of an increase in powerplant weight, since 
the structural weight changed very little. The lift-drag ratio (tables II and m) decreased 
from about 17 to 15, which corresponded to an 11-percent drag increase. Although the 
total drag coefficient decreased 25 percent, mostly because of lower drag due to lift, the 
product of dynamic pressure and wing area 1 qS w increased 50 percent, mostly due to 
Mach number, and thus resulted in the previously quoted drag increase. Similar rea- 
soning should apply to the 45 000-foot (13 720-m) case. 

Calculations over a broad range of off- optimum wing loadings at Mach 0. 8 from 60 

o 

to 119 pounds per square foot (287. 2 to 578. 5 kg/m ) showed that the payload was, in- 

2 

deed, maximum for the indicated optimum wing loading of 85.2 (408 kg/m ) (cruise lift 


25 






coefficient of 0.401). For the specified fixed conditions (Wg, Mq, H, C^), the weights 
of the powerplant and wing are related through aspect ratio but not in a simple or direct 
manner. Thus, in effect, aspect ratio was varied by the optimizing procedure to find 
the value that gave the best payload (the other independent variables were free to seek 
new values if profitable). 

These results should be considered preliminary since additional considerations of a 
more detailed nature, for example, wing flutter, may alter the trade-offs or add other 
boundaries. 

The cruise lift-drag ratios for the spectrum of design cases were as high as 18. 3, 
but decreased to between 16 and 14 for the highest Mach number considered for each al- 
titude. The Mach number for maximum payload, figure 5(c), did not correlate with 
maximum lift- drag ratio or wing loading. 

The total stream tube areas are presented in figure 5(c) primarily to illustrate the 
variation of design-point engine airflow requirements. Furthermore, if a velocity ratio 
of 1. 0 is selected between the stream and the blading (Vj/Vq), a fan diameter can be es- 
timated. For a six engine airplane, the largest area corresponds to a fan diameter of 
about 13 feet (3. 96 m) and the smallest is about 7. 4 feet (2.25 m). This estimate of fan 
diameter only affects the calculated results by way of the nacelle weight and drag which 
are a very small fraction of the total drag and structural weight (appendixes D and E). 

Some other very general trends can be observed from the data of tables II and III. 
The installed power varied from about 175 to 384 megawatts, of which from 5 to as much 
as 10. 4 percent was needed for helium pump power. For the engine cycle parameters 
the overall compressor pressure ratio was between 14 and 21, and it generally reached 
a peak at some intermediate value of Mach number. Fan pressure ratio increased with 
flight Mach number and altitude, and was within the range 1. 2 to 1.6. Bypass ratio var- 
ied between 8. 4 and 2. 9, and decreased with increasing Mach number as cycle efficiency 
increased. Turbine inlet temperatures were in the range 1800° to 2100° R (1000 to 
1166 K), and they generally were higher at the highest altitude. 

Since airplanes designed for altitude extremes such as 20 000 or 45 000 feet (6096 
or 13 720 m) may be operationally limited, an altitude of 36 089 feet (11 000 m) was se- 
lected for the remainder of this study. For example, the usual "over-weather” crite- 
rion requires an altitude capability of at least 30 000 feet (9140 m), and flight at an alti- 
tude of 45 000 feet (13 720 m) is no advantage for logistic missions. A Mach number of 
0. 8 was also selected as a base point, when necessary, since most advanced subsonic 
cargo planes fly at or near this speed. 

Effect of gross weight. - Payload fraction increases with gross weight at a decreas- 
ing rate for either Mach 0. 6 or 0. 8, as shown in figure 6. This occurs because the de- 
crease in powerplant weight fraction with gross weight has leveled out and the, structural 
fraction increased a few percent of gross weight at the higher gross weights. Delivery 
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Gross weight, Wg, kg 

Figure 6. - Variation of weight fraction and delivery rate with gross weight for 
the optimized system. Design altitude, 36 089 feet (11 000 m) : shield, metal- 
water with packing factor of 4/3. 


rate also increases monotonically with gross weight. The magnitude of the payload for 
the 1/2 -million- pound (22. 7xl0^-kg) airplane is one-eighth of that for the 1-million- 
pound (45. 4x10 -kg) aircraft. Again, these trends are representative of the parametric 
families and technology level used. The transition from preliminary to detailed or lay- 
out design may reveal exceptions or alternates to the rules employed; for example, new 
materials (such as composites) and structural techniques may favor large aircraft. The 
variation of the weights of components of the powerplant and structural groups with air- 
craft gross weight is presented in figure 7 for the base-point conditions for the optimized 
system at Mach 0. 8 and 36 089 feet (11 000 m). Of particular interest as gross weight 
increases is the counteracting effect between the favorable decrease in shield weight (in- 
cluding duct penetration) fraction (W Cin /W„ = 0. 39 to 0. 17) and the increase in the 
chemical fuel, engines, heat exchangers, and lines. The sum of these latter items in- 
creased from Wgjj M /Wg = 0. 18 to 0.22 and thus represents an important weight group. 
For the lower gross weight, the shield weight fraction is proportionally greater (than for 






.2 .3 .4 .5 . 6x10*’ 

Gross weight, Wg, kg 

Figure 7. - Detailed weight breakdown as function of gross weight. Design 
flight Mach number, 0.8; altitude, 36 089 feet 111 000 m); shield, metal- 
water with packing factor of 4/3. 


higher W g 's) not only due to the lower installed power (eq. (1)), but also due to the de- 
creased dose distance (eq. (2)) and the constant dose level criterion. 

The relative percentages of the weights of various items in the structural group 
(such as wing, fuselage, and tail) were also relatively constant over the range of gross 
weights (fig. 7) except for those items specifically programmed as functions of gross 
weight (such as the fixed equipment and systems, and landing gear). 

The reason for the chemical fuel fraction increasing with gross weight is related to 
the general decline of engine efficiency and lift-drag ratio (as increases). Again, 
the trade-off is between installed power and wing area as given in the previous section. 
However, because the power changes with gross weight the program seeks a balance be- 
tween the shield weight as affected by reactor flow rate w^ e /A and helium temperature 
rise (different relative core size) which, in turn, varies the helium pressure drop and 
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hence pump power requirements. Extracting pump power from the engine cycle changes 
the overall engine efficiency and, therefore, alters the chemical fuel allowance. The 
engine efficiency during chemical operation, which is without pump power extraction, is 
higher than during nuclear operation, but follows the trend for the nuclear operation 
since the cycle parameters are optimized for the nuclear cruise condition. 

The lift-drag ratio varies (slightly) as a result of obtaining the best balance between 
wing weight (size and geometry) and powerplant weight for maximum payload. 

Effect of engine type . - The optimum propulsive system defined in the preceding re- 
sults, wherein bypass ratio was one of the twelve independent variables, was in all cases 
a turbofan. A comparison of the performance and nominal values of certain cycle pa- 
rameters for the optimized system using turbofan (BPR > 0) or turbojet (BPR = 0) en- 
gines is shown in figures 8 and 9. The advantage for turbofan propulsion in terms of 
payload fraction varied from about 16 percent improvement at Mach 0. 8 to almost 
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Figure 8, - Comparison of turbojet and turbofan propulsion for optimized 
system. Gross weight, 1 million pounds (453 600 kg); altitude, 36 089 
feet; (11 000 m) shield, metal-water with packing factor of 413. 
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Design flight Mach number, Mq 


Figure 9. - Comparison of turbojet and turbofan cycle parameters for 
optimized system. Gross weight, 1 million pounds (453 600 kg); 
altitude, 36 089 feet (11 000 m); shield, metal -water with packing 
factor of 4/3. 


108 percent at Mach 0. 4. This is primarily related to the inherent higher overall effi- 
ciency of the turbofan cycle at lower Mach numbers. The fan cycle was nearly twice as 
efficient as the jet cycle at this low speed condition and was 5 to 8 efficiency counts 
higher over the range of Mach numbers (fig. 9). The jet cycle attempts to counter this 
effect of lower overall efficiency by seeking higher cycle pressure ratios but is tempered 
by the effect of increasing engine weight. The fan cycle favors higher turbine inlet tem- 
peratures primarily because of the need for turbine power to drive the fan. 


Effect of Off -Optimum Component Performance 


Two general types of component performance variations are classified herein. The 
first type considers large changes such as advances in technology or major alterations 
of initial goals. The second type referred to as sensitivity factors is for smaller incre- 
ments of the component performance such as may occur during detailed design. 

Effect of overall compressor pressure ratio . - The importance of overall compres- 
sor pressure ratio to payload and the accompanying values of various other parameters 

4 

are shown in figure 10 for the 1- million- pound (45. 4x10 -kg) airplane at Mach 0. 6 and 
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0. 8. The optimum overall pressure ratios which varied with flight design point have 
been presented in tables n and m. 

Overall compressor pressure ratio affects payload primarily through the weight of 
the installed powerplant. The installed power varies inversely with the overall effi- 
ciency as shown by 

V 0^-”a 

V (f 


In appendix A, the engine weight is given as a function of overall pressure ratio, by- 
pass ratio, turbine temperature, and thrust. This, together with equation (8), illus- 
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Turbine inlet temperature, T 4 , K 

Figure 11. - Effect of turbine inlet temperature on weight distribution. Gross weight, 1 million 
pounds (453 600 kg); design Mach number, 0.8; altitude, 36 089 feet (11 000 m); shield, metal-water 
with packing factor of 4/3. 
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TABLE IV. - EFFECT OF HIGHER THAN OPTIMUM TURBINE 


INLET TEMPERATURE FOR BASE-POINT CONFIGURATION 



Turbine inlet temperature, T^ 

Percent 

change 

a 1890° R (1049 K) 

2250° R (1249 K) 

Weight fraction 




Heat exchanger 

0.048 

0. 0922 

92 

Engine 

.063 

. 050 

-20 

Ducts 

.015 

. 020 

33 

Payload 

. 192 

. 145 

-25 


a This temperature corresponds to maximum W^/W^. 


trates the trade-off between shield weight (function of eq. (1)) and engine weight for 
maximum payload as overall compressor pressure ratio is varied. Values of the other 
independent variables for the engine that corresponded to maximum payload for a pre- 
scribed overall compressor pressure ratio are shown in figure 10(b). 

The symbol representing the ” computer selected” independent variable values for 
best payload do not always lie on the faired lines of parametric values. This small dis- 
crepancy represents a band of alternate selection within which the independent variables 
may be mutually traded without altering the payload (within its convergence tolerance). 

Effect of turbine temperature . - The effect on payload fraction and some engine 
cycle parameters of varying the turbine inlet temperature over the range 1800° to 
2250° R (1000 to 1249 K) is shown in figures 11 and 12. The nominal maximum reactor 
wall temperature was 2560° R (1422 K) and the temperature of the helium leaving the 
reactor was not permitted to exceed 2360° R (1311 K). Because of the freedom to re- 
optimize other components, the payload fraction varied only 7 percent over the range of 
turbine temperature of 1800° to 2100° R (1000 to 1166 K), but was 25 percent below the 
maximum at 2250° R (1249 K). The peak payload value occurred at 1890° R (1049 K). 

At the higher temperatures, the heat exchanger weight increased more than the engine 
weight decreased as shown in the powerplant weight breakdown and in table IV. For the 
lower turbine temperatures, the cycle parameters are generally reduced due to the dif- 
ficulty of maintaining a positive work out of the turbine. 

The small payload penalty for the lower turbine inlet temperatures suggests that 
lower reactor wall and helium temperatures, having simpler materials problems, may 
give comparable results. 

Effect of heat exchanger pressure drop and lifetime . - The heat exchanger is the in- 
terface between the helium and air sides of the system. Once the heat exchanger is pro- 
portioned to transfer the required energy to the air, the resulting air-side pressure 
drop, although the magnitude is important to the cycle performance, assumes an almost 











constant value of 10 percent (table II) (includes burner pressure loss) and is apparently 
not strongly affected by the feedback effect of engine weight on the helium-side weight. 

An air-side pressure drop of this magnitude causes about a 10-percent reduction in over- 
all efficiency n „ . 

Extending the calculated lifetime of the heat exchanger from 1000 to 10 000 hours 
had a very small effect on payload. Heat exchanger weight increased 9250 pounds 
(4195 kg), or 19 percent, and payload decreased about 9800 pounds (4443 kg), or 
5. 1 percent, for the base-point configuration (Mq = 0. 8, altitude = 36 089 ft (10 990 m), 
W g /10 6 = 1. 0 lb (0. 4536 kg), cp = 4/3). 

Effect of system pressures . - The effect of changing the helium pressure from the 
value of 1750 psia (12. 06x10)6 N/m^) thus far considered is shown in figure 13. No max- 
imum was found within the range from 1600 to 2000 psia (11. 03 to 13. 8x10^ N/m^), and 
the increase in payload over the range was only 1. 5 percent. 
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Figure 13. - Effect of helium system pressure. Gross weight, lxlO 6 
pounds (453 600 kg); design Mach number, 0.8; altitude, 36 089 
feet (11 000 m) ; shield, metal -water with packing factor of 4/?. 


As shown in the figure, the higher helium pressures required slightly less installed 
power, had less pressure drop ratio, and hence less shaft power extraction from the en- 
gines. This allowed a higher efficiency for the engine cycle and hence less chemical fuel 
was required which accounted for about one -third of the small improvement in payload. 
The remainder of the improvement was primarily concentrated in the smaller shield, heat 
exchangers, pumps, and reactor weights associated with the smaller sizes at the higher 
system pressures. The small increases in duct and engine weights caused by increasing 
system pressure were easily counter-balanced by the previously mentioned factors. 

The variation with Mach number of the total system helium pressure drop ratios and 
the distribution by components is presented in figure 14 for the nominal helium pressure 





Figu re 14, - Va riation of component pressu re drop ratios with design 
flight Mach number. Gross weight, lx# pounds (453 600 kg) ; 
altitude, 36 089 feet (11 000 m); shield, metal -water with packing 
factor of 4/3. 

o 

of 1750 psia (12. 06 N/m ). The accompanying variation in pump and installed power is 
given in tables II and in. (It should be noted that the component pressure drop ratios 
are not additive in the form presented since each is based on the entrance pressure for 
that component. ) The increase in system pressure drop ratio from 0. 075 to 0. 124 over 
the Mach number range of 0. 4 to 0. 8 is primarily due to the increases for the reactor 
and shield components. This occurs as the installed power increases with Mach number 
(see tables n and HI). There is a trade-off between the savings in powerplant weight for 
smaller reactor radius (shield weight effect) and the increased pumping requirements. 
The pressure drop ratios shown represent the best combination that maximizes payload 
and not necessarily the minimum pressure drop for any one component. This is also an 
example of the reason that performance-weight relations are needed for the major com- 
ponents. 

Effect of specific engine weight. - The effect of improving the specific engine weight 
is shown in figure 15 for the 1-million-pound (45. 4x10 -kg) airplane. The system per- 
formance has been reoptimized. A significant improvement of about 28 000 pounds 
(12 690 kg) of payload (2. 8 percent of gross weight) was obtained for either the turbofan 
or turbojet at Mach 0. 6 or 0. 8 for a reduction in engine weight of 40 percent (engine 
weight is about 6 percent of the gross weight). 

Effect of fuselage length and reactor shield fairing . - As discussed in a previous 
section, some of the high altitude and/or Mach number design cases (fig. 5(a)) had such 




Ratio of improved to base-point 
engine weight 

Figure 15. - Effect of engine technology on perform- 
ance of the optimized system. Gross weight, lxlO 6 
pounds (453 600 kg); altitude, 36 089 feet (11 000 m); 
shield, metal -water with packing factor of 4/3. 


low payload fractions that the fuselage size exceeded the need. Therefore, the effects of 
placing the reactor within the cargo bay (thus restricting the maximum length of cargo 
items and drive-through capability) and also of shortening the fuselage were investigated. 
Results are presented in figure 16. 

Over the limited range of conditions investigated, the payload gains realized by re- 
ducing the size of the reactor-shield fairing were generally less than 4000 pounds 
(1814 kg) for the 1- million- pound (45. 4xl0^-kg) airplane. This is a trade-off between a 
small reduction in installed power due to slightly improved lift-drag ratio (less surface 



Figure 16. - Effect of fuselage size changes on performance of 
optimized system. Gross weight, 1x10° pounds (453 600 kg); 
shield, metal-water with a packing factor of 4/3. 
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(a) Payload and nuclear side parameters. - (b) Engine cycle parameters. 

■Figure 17. - Effect of reactor helium weight flow per unit area on system performance. Gross weight, 1x10*’ pounds (453 600 kg); design Mach number, 0.8; 
altitude, 36 089 feet (11 000 m); shield, metal-water with packing factor of 4/3. 









area) and some increase in fuselage weight as the section moment of inertia is decreased 
(eq. (E16)). 

Decreasing the fuselage length to width ratio from 12 to 9 resulted in a payload in- 
crease of between 20 000 and 30 000 pounds (9070 and 13 600 kg). The fuselage profile 
drag was reduced about 25 percent (or about 7-percent reduction in airplane profile drag) 
thus increasing the lift-drag ratio about three-tenths of a unit and concommitantly de- 
creasing the installed reactor power. As the fuselage is shortened, the dose distance 
decreases from 130 to 97. 5 feet (39. 6 to 29. 72 m) and increases the shield weight frac- 
tion, according to equation (2), about 10 percent, or 20 000 pounds (9070 kg). This 
shield weight rise is balanced, in general, against a fuselage weight decline of 40 000 
pounds (18 140 kg) to arrive at the cited advantage. 

Effect of reactor helium flow rate per unit area . - A typical performance map for a 
family of reactor designs has been shown in figure 2. The power level and helium weight 
flow per unit area (of coolant passage) are used to establish both a core radius and a 
correlation factor used in the shield weight calculations (see eq. (1)). The choice of 
using a higher reactor flow rate per unit area to obtain a smaller and lighter shield must 
be balanced against the consequences of an increased pressure drop. 

As shown in figure 17(a), payload fraction increased with reactor helium flow rate 
per unit area only about 5 percent or 1 percent of gross weight over the range considered 
(w He /A of 120 to 180 (586 to 879 kg/(m 2 )(sec)); AP/P He reactor 3 to 7^ percent) at 
Mach 0. 8 and 36 089 feet (11 000 m) in altitude. Notice that the shield weight fraction 
decreases about 2 percent of gross weight as Wjjg/ A and installed power increase be- 
cause of the smaller core size. Installed power increases for the following reason: In- 
creasing the reactor flow rate parameter w^/A causes an increase in the core pres- 
sure drop (fig. 17(a)), and hence the total pressure drop, which increases the pump 
power required, and because the pump is shaft powered, engine efficiency decreases 
(fig. 17(b)). This decrease in engine efficiency increased the chemical fuel allowance 
about 1 percent of gross weight. Other components of the powerplant group, in addition 
to the shield and chemical fuel underwent readjustment, with the net result of very little 
change in this weight fraction. However, the structural weight fraction was reduced 
slightly more than 1 percent of gross weight primarily due to a lighter wing. Thus, the 
final result was dependent on interactions among many components, primarily following 
the network previously established such as in the section Effect of gross weight (p. 26). 

Shield technology . - Thus far, all the results presented have been for the metal- 
water "base-point" shield with a packing factor (or shield weight multiple) of 4/3. The 
difference between the metal-water and hydride-water shields and the meaning of the 
term packing factor are discussed in the section Shield (p. 12). The variation of payload 
fraction with packing factor for both shield designs is shown in figure 18. Packing factor 
and shield technology, both representing substantial adjustments in shield weight, were 




Figure 18. - Variation of payload weight frac- 
tion with packing factor and shield tech- 
nology. Gross weight, 1x10^ pounds 
(453 600 kg); design Mach number, 0.8; 
altitude, 36 089 feet (11 000 m). 
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Figure 19. - Influence of shield technology on major weight fraction as func- 
tion of gross weight. Design Mach number, 0.8; altitude, 36 089 feet 
(11000 m); shield packing factor, 4/3. 


very important. In fact, the combination of advanced shield materials (hydrides) and a 
very small packing factor («1. 0) could increase the payload fraction by a factor of 
50 percent compared with the previous result for metal-water shields having a packing 
factor of 4/3. The improvement in payload fraction due to shield technology is shown in 
figure 19 to be very beneficial over the range of gross weights from 1/2 to 1-^ million 
pounds (226 700 to 566 500 kg) and more important at the lower gross weights where 
powerplant weight components, the shield in particular, account for the general decrease 
in payload fraction. 

With advanced shield technology, the payload fraction of 16 percent obtained for the . 
1/2 -million-pound (226 700-kg) airplane is about three times that for the metal-water 
shield with a packing factor of 4/3; however, at Wg = 10° pounds (453 600 kg) the ad- 
vantage decreases to 1.3. 

Component sensitivity factors . - Two concepts of sensitivity factors are possible 
depending on whether or not the design is considered frozen. If the system design is 
frozen and the performance of a certain component must be revised, the design (oper- 
ating points) of the remaining components is not altered, and the term ’’fixed" is applied. 
On the other hand, if the revised performance of a component can be offset by redesign- 
ing or reoptimizing the other components, the term "reoptimized" or "rubber" is used. 

Examples of each type of sensitivity factor are given in table V for each of the inde-- 
pendent variables (listed in Optimization Procedure (p. 19) and appendix F). Where 
possible, a ±10 percent change was used. The data are for the base-point case of 
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Variable 


Percent change 
in variation 1 


Fixed 


Rubber 


Percent decrease in payload 


P9/P1 


10 

-10 


5 . 17 
3.38 


0 . 52 
. 67 


V P 2 


10 

-10 


2.2 

4.4 


0.26 
. 39 


P 4a/ P 3 


2 . 85 
- 2 . 85 


0.27 
. 69 


0.29 

1.02 


BPR 


10 

-10 


0.22 

.314 


0.23 
. 306 


AR 


10 

-10 


0 . 32 
. 81 


0.28 

.79 


10 

-10 


0.8 

. 581 


0 . 63 
. 579 


10 

-5 


12 . 54 
3 . 8 


5 . 77 
1 . 53 


AP. 


HX 


10 

-10 


0 

0 


AP. 


ducts 


10 

-10 


W He/ A 


-10 


1. 12 


1.43 


l HX, out 


10 

-10 


8 . 0 

7 . 6 


4.7 

2 . 6 


L He 


-10 


11 . 6 


7.3 


H 


10 

-10 


13 . 75 
b 3 . 4 


13.9 
b 3 . 88 


M r 


10 

-10 


21 . 94 
b 9 . 4 


18 . 97 
3 15 . 5 


d Any change less than tolerance of convergence is 
considered zero. 
b Gain. 
































































M n = 0. 8 , altitude of 36 089 feet (11 000 m), metal-water shield, packing factor of 4/3, 

U a 4 

and Wg = 10 pounds (45. 4x10 kg). Since the system was optimized, all of the sensi- 
tivity factors are negative or zero if the result was less than the computing tolerance. 

The independent variables that were most sensitive to change were, in order of rel- 
ative importance, turbine inlet temperature, helium temperature out of reactor, helium 
temperature out of heat exchanger, and overall and compressor pressure ratios. Al- 
lowing reoptimization decreased the sensitivity factors; for instance, those involving 
temperature were reduced by about 50 percent. For a few independent variables (neg- 
ative P 4 /P 3 values and w 

or " rubber” case and was the result of the selected convergence tolerance on payload. 


He /A) a greater sensitivity is indicated for the reoptimized 


Comparison With Conventional Aircraft 

The estimated performance of the nuclear aircraft is compared in figure 20 with the 
payload- range values predicted for the Air Force C-5A airplane. (The C-5A is the larg- 
est logistics aircraft under construction at this time. ) Data are presented for both one- 
way and round trips, no refueling missions. Gross weights of 1.25, 1. 0, and 0.73 mil- 
lion pounds (56. 7xl0 4 , 45. 4xl0 4 , and 33. lxlO 4 kg) are shown for the nuclear plane, the 
latter value being approximately equal to that for the C-5A. For round-trip flights, the 
chemical fuel for an additional takeoff and landing at the midway point has been deducted 
from the payload for the nuclear planes. The flight Mach number of 0. 8 for the nuclear 
plane is approximately equal to that for the C-5A. 

For equal gross weights, the nuclear plane has superior performance for ranges 
greater than about 5900 nautical miles (~llxl0° m) for one-way missions or about 2700 

miles (~5xl0 m) for no- refueling round trips with payload delivery at the midway point. 

4 

A 1- million-pound (45. 4x10 kg) nuclear airplane would exceed the payload of the C-5A 
at ranges beyond 3700 nautical miles ( 6 . 9x10° m) for one-way trips or beyond about 2000 

nautical miles (3.7x10 m) for round trips. Payload estimates for the l4- million- pound 

4 ^ 

(56. 7x10 -kg) nuclear plane exceed those for the C-5A at all ranges of interest for either 

one-way or round trips. The nuclear airplane's virtually unrestricted range (or, more 
correctly, endurance) is only limited by human capabilities and equipment reliability. 

Flight and ground handling considerations are not comparable for the two systems. 
However, the use of a unit shield and low dose rates averts some of the operational 
problems encountered in prior nuclear airplane concepts. No penalties have been as- 
sessed for crash safety provisions such as energy absorption systems for the nuclear 
airplane. 

Technology advances in aerodynamics, propulsion, and structures would be mutually 
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Figure 20. - Comparison of nuclear and conventional fossil fueled aircraft on 
the basis of payload and range. Nuclear aircraft's design flight Mach number, 

0.8; altitude, 36 089 feet (11 000 m); shield, metal -water with a packing factor 
of 4/3. 

beneficial for either the nuclear or conventional planes, whereas advances in nuclear 
technology would, of course, only benefit the nuclear concept. 

SUMMARY OF RESULTS 

A preliminary design-point investigation of a subsonic logistic aircraft powered 
during cruise by a helium-cooled unit-shielded thermal nuclear reactor has been con- 
ducted. A range of flight design points, engine cycle parameters for both turbojets and 
turbofans, reactor coolant flow rates, and pressure have been considered. Variations 
in gross weight and shield and engine technology were studied, as well as the interrela- 
tions and sensitivity factors of the major components. The maximum helium tempera- 
ture was 2360° R (1311 K) and the nominal maximum reactor wall temperature was 
2560° R (1422 K). For the base-point or standard calculation, the helium pressure was 
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1750 pounds per square inch (12. 06x10 N/m ), and a metal-water shield with a fabri- 
cated over ideal weight ratio of 4/3 (called packing factor) was used. The radiation dose 
level was constant at the very low value of 2. 5 millirem per hour at the flight deck, and 
for 1000 hours of flight duty corresponds to only one-half the man-year recommendation 
of the Federal Radiation Council. Payload capacity was selected as the figure of merit 
to be maximized. The following results were obtained: 

4 

1. For the 1-million-pound (45. 4x10 -kg) gross weight state-of-the-art configura- 
tion, the payload fraction was 19 percent for a design altitude of 36 089 feet (11 000 m) 
at Mach 0. 8. The effect of increasing design altitude and Mach number was to decrease 
the payload. At lower design altitudes and Mach numbers (20 000 ft (6096 m) and 0. 5) 
payload fractions increased to about 26 percent; however, such ceiling limited airplanes 
may be of limited operational value. For the Mach 0. 8, 36 089-foot (11 000- m) design 
point, the turbofan bypass ratio was 3. 6, the overall pressure ratio was 18. 4, the tur- 
bine inlet temperature was 1890° R (1049 K), and the pressure drop between the com- 
pressor and turbine due to the heat exchanger was 12 percent. On the nuclear side of 
the system, the installed power was 309 megawatts (including 9 percent for pump power), 
the helium pressure was 1750 psi (12. 06x10 N/m ), and the pressure drop was 12 per- 
cent. The airplane had an optimum wing loading of 85 pounds per square foot 

o 

(415 kg/m ) and a structural weight fraction of 36. 7 percent including equipment and 
systems. 

2. Since shield weight represented about 20 percent of the gross weight for the 1- 

4 

million-pound (45. 4x10 -kg) airplane, shield technology had an important effect. At a 
design Mach number of 0. 8 and an altitude of 36 089 feet (11 000 m), the payload fraction 
could be increased 30 to 50 percent by using a hydride-water shield material rather than 
the metal- water type used as the standard case. 

3. Payload fraction increased with gross weight at a declining rate up to the largest 

considered. The value for the 1/2 -million- pound (226 700-kg) airplane was 5 percent 

1 4 

and increased to 19 and 22 percent, respectively, for 1- and 1^- million- pound (45.4x10 - 
and 56. 7xl0 4 ~kg) airplanes designed for Mach 0. 8 and 36 089 feet (11 000 m). Although 
the shield was still the predominate weight fraction component and decreased as gross 
weight increased, the weight fractions of the structure and a group composed of the en- 
gines, heat exchangers, and chemical fuel for nonnuclear operation increased. 

4. The optimum turbine inlet temperature for maximum payload was within the 
range of 1800° to 2100° R (1000 to 1249 K), depending on the flight design point. For 
Mq = 0. 8 and 36 089 feet (10 990 m), a turbine inlet temperature of 1890° R (1049 K) 
maximized payload which varied only 7 percent for turbine temperatures between 1800° 
and 2100° R (1000 to 1249 K) because the other components were allowed to reoptimize. 
At these higher than optimum turbine temperatures, heat exchanger weight was increas- 
ing more rapidly than the engine weight decreased. 



5. Varying the system pressure from 1600 to 2000 psia (11. 03x10® to 13.78x10® 

2 

N/m ) increased the payload fraction only about 1. 5 percent, mainly due to reductions in 
installed and pump power. 

6. The payload fraction obtained with the optimized system using the turbofan cycle 
was 16 percent greater than for the turbojet cycle at Mach 0. 8, and 108 percent greater 
at Mach 0. 4 for a design altitude of 36 089 feet (11 000 m). The overall efficiency for 
the turbofan cycle was 5 or 8 efficiency counts higher, while the optimum overall com- 
pressor pressure ratio was higher for the turbojet. The optimum turbine temperature 
for the turbofan was generally 100° to 200° F (55. 5 to 111 K) higher than for the turbojet. 

7. Sensitivity effects due to perturbing the independent variables were more pro- 
nounced when considered individually (fixed design) than when other independent vari- 
ables were allowed to reoptimize. In general, the off-optimum temperature effects 
were more detrimental. 

5 

8. Comparison with the C-5A class of airplanes at equal gross weights (~7. 3x10 lb; 
331 000 kg) indicates superior payloads for the nuclear plane for ranges beyond 5900 
nautical miles (~llxl0° m) for one-way trips or 2700 nautical miles (~5xl0° m) for 
round trips without refueling. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, April 1, 1968, 

789-50-01-01-22. 
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APPENDIX A 


TURBOFAN-JET WEIGHT AND COMPONENT ASSUMPTIONS 


The engine weight equation was taken to be of the form 


f W, 


w 


eng 


eng' 


SLS 


r SLS 


where the specific weight at sea level static (W. _/F) was assumed to be a function 

enr SLS 

of overall pressure ratio, bypass ratio, and technology level. The sea level static 
thrust Fg^g is computed from the cruise thrust and factors that account for design al- 
titude, Mach number, and turbine inlet temperature. In particular, 


W \ /P„\ W 

6ng ' = f i ( — Vo ( BPR ) f 3 (T 4 ) — 


r SLS 


where 




0. 834 exp 


0. 722 exp 


v. 


^no. o ) 

55. 5 


(2 -spool engines) 


(1-spool engines) 


f 2 (BPR) = 0. 97 - 0. 001 (BPR) + 0. 0014 (BPR; 
f 3 (T 4 ) = [l + 5. 5xl0" 4 (T 4 - 2060)] _1 (T 4 in °R) 


W 


= reference specific weight (function of technology 
F level), assumed to be 0.213 at P^/P^ = 20, 

BPR = 5, and T 4 = 2060° R 
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W'* = ^ at 36 000 ft) 


F SLS = Vf 
F des a 


1 + M r 


% - t *) 2 
a - *7 


The last relation is not valid if BPR > 8 or for altitudes less than 20 000 feet if 
BPR <3. 

Although F sLS^ F des is also a functi on of BPR, its dependence on BPR is quite 
small at the high design altitudes (e. g. , 36 000 ft) of major interest here. It is also a 
good approximation at lower design altitudes when BPR ~ 4. These f factors and the 
F SLS^ F des functl0n were selected after compiling data from various engine designs and 
prototype studies and are displayed in figure 21. Table VI gives the component efficien- 
cies assumed for this study. 


TABLE VI. - COMPONENT EFFICIENCIES AND 


PRESSURE RATIOS 


Component 

Assumed efficiency, 
percent 

Fan 

0. 88 (adiabatic) 

Compressor 

. 86 (adiabatic) 

Turbine (fan) 

. 90 (adiabatic) 

Turbine (compressor) 

. 90 (adiabatic) 

Fan duct thrust coefficient 

. 985 

Primary nozzle thrust coefficient 

.985 

Pressure ratio in primary nozzle 

. 98 

Inlet pressure recovery 

1.0 

Pressure ratio in duct nozzle 

.97 

Pressure ratio between turbine inlet 

. 95 

and compressor outlet exclusive 


of heat exchanger 

_ 
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APPENDIX B 


TAYLOR EXPANSION FOR REACTOR PERFORMANCE MAP AND 
EQUATIONS FOR SHIELD WEIGHTS AND RADIUS 

Taylor expansions were used to curve fit results from the reactor analysis to relate 
the reactor power, pressure drop, helium temperature, radius, and helium weight flow 
per unit area. The following terms are used to simplify the expressions developed: 


AT- = T. - 

T„ 

in in 

in, o 

AT = T 

T 

ex ex 

ex, o 

AP. = P. - 

P. „ 

in in 

in, o 


A& = 9 


AW _ W 
A A 



where the nominal maximum reactor wall temperature is constant at 2560° R, T in is 
the helium pump helium exit temperature (°R), T ov the reactor helium exit temperature 
(°R), P in the reactor inlet helium pressure (psia), 9 the reactor power (MW), and 
W/A the reactor helium flowrate per unit free area (lb/sec/ft^). The subscript o de- 
notes the reference point values of 1260° R, 2260° R, 1500 psia, 250 megawatts, and 
120 pounds per second per square foot, respectively (equivalent to w He /A in main body 
of report). 

The expansion for pressure drop can be written in a simplified manner by defining 


9.2696x10 ~ 3 (~^ 


g x = exp 

g 2 = exp(-2.215xlO- 3 P ino ) 


Then, 


AP 


rea 


0. 0145 g 1 - 0. 02 - 1. 027x10" 5 AT. + (2. 24xl0" 7 T - 4. 427xlO" 4 )AT 

J- 111 vAj U CA 


-4 „ AW 


+ 1. 12xlO" 7 (AT ex ) 2 + 1. 345X10" 4 g 1 
+ 6.25X10" 7 g 1 - 1. 0034X10 -3 g 2 AP in 

+ 1. HlXlO" 6 g 2 (AP in ) 2 - 8. 2033x10" 10 g 2 (AP in ) 3 


1. 3889 


The reactor radius was calculated from the Taylor expansion: 


R„. = 1.42206x10' 

led. 


>/w\ 2 _ 

\ A/ 


6. 9255x10" 3 (— ) + 2. 5973 + 3. 1593X10" 3 A& 


2. 84412x10 


-5 


- 6. 9255x10' 


+ 1. 42206xl0~ 5 ^— J 


The T- and T terms do not appear in this expression because their effect on 
in ex 

R was negligible as determined from detailed calculations over the range of tempera- 
tures used herein. 

Reactor weight was calculated from 

= 24. 72 9 + 682 

icd 


which was then doubled to account for moderator (water) to air heat exchanger and after- 
heat removal system weights. 

A curve fit to f, which is in reality a measure of the reactor core density, was then 
found: 


R„ oq (in cm) = 30. 48 R (in ft) 

red. 

R 1 = -10" 5 ^ 2 + 0. 089^ + 44. 9750 
P ?2 = -3. 5xl0" 5 ^ 2 + 0. 097 + 48. 6375 
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Kg = -4. 5x10' V + 106^ + 51. 5125 
C 1 =H |- R2 r e a< incm ) 


C 2 ~ **2 " R 1 



C 4 “ R 3 " *2 


X = 


C 3 ~ C 1 
C 3 C 2 ' C 1 C 4 


Y = 


i-c 2 b 


Z = 1 - YR“ - XRj 


Then, 


f = YR^o (in cm) + XR rea {in cm) + Z 


rea 


If f < -4, set f = -4. 

For base-point shields (metal-water), 

JO. 122 

R„ = (113. 5 + 2. 5 f) — 

s 30. 48 

W = (17 900 + llOOf)^ 0 ’ 407 

o 

For advanced shields (hydride- water), 


R = (85 + 2.0 f)£ 


0. 134 


30.48 


W e = (12 100 + 1000f)^ 
s 


0. 419 



APPENDIX C 

CHEMICAL FUEL ALLOWANCE 

The amount of fuel required for the chemical flight phases is primarily dependent on 
the engine specific fuel consumption (sfc), the basic design altitude and Mach number, 
and the emergency range (or duration). To compute the fuel allowance, the total fuel re- 
quirement was broken down into three parts: ascent, cruise, and descent. Thus, 

W F = W F , T c +W F, c + W F,d ( C1 ) 


Ascent Fuel Wp jq 

The climb fuel requirement can be approximated by energy balance relations. The 
change in energy per pound of aircraft weight (specific energy) may be written (ref. 10, 
eq. 2) 


Wde = 7? oa HV dW p> TC - DV q dt (C2) 

where W is the aircraft weight, e its specific energy, r) Qa the engine overall effi- 
ciency, HV the fuel heating value, D the airplane drag, V Q the flight velocity, and t 
the time. Since dH = Vq sin a dt, where a is the path angle, equation (C2) may be re- 
arranged and formally integrated to give 



T, oa HV dW 


F, TC ~ 



(C3) 


This equation is still exact; that is, it accounts for the energy change due to change 
in gross weight as fuel is burned. Defining rf to be the average efficiency and W to 
be the average weight during climb, equation (C3) can be solved for the climb fuel: 


W 


F, TC 



(C4) 
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For a first-order approximation, W can be taken to be the same as the takeoff 
weight. Also, 77 oa can be considered to be independent of altitude but varying from zero 
at V = 0 to cruise efficiency Tj Qa c at V = Vq. Thus, for a change in altitude only, 

77 would be taken to be But for an increase in velocity from zero to V n , rf~Z 

would be taken to be 1/2 « . And if during climb both velocity and altitude change, 

then 7 j oa could be calculated as the sum of two weighted terms - one term associated 
with the change in kinetic energy KE and the other potential energy PE. For simplic- 
ity, the function chosen was 


^oa 


/ 1 AKE + 
\2 AE 


1 . 0 


APE\ 

ae r° a > c 


(C5) 


where AE is the change in total aircraft energy during climb (i. e. , AE = W Ae). 
The integral of D/sin a was approximated by rewriting it as: 



C D^ pV 0 S w C D S 
dH = D w 


sin a 


2 sin a 


PV 2 dH 


(C6) 


where again bars denote flight path averages. The integral of pVq can be approximated 
by assuming that the climb path is a straight line on a M-H diagram between the points 
(0. 2 5, 0) and (M d , H d ). Then, 



M = aH + b 


where 


M d - 0.25 
H d 

b = 0.25 


Also 


Vq = M 2 yRT (y = 1.4) 
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where 

(H d < 36 089) 

(H d > 36 089) 

The second term of equation (C7) is not required for H > 36 089. All that remains is to 
evaluate C^/sino. Letting C D d be the cruise design value of gives 
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(C8) 



sin a 


'D, d 


'D 


'D, d 


28. 65 C 


sin a 


D, d 


where the constant 28. 65 was chosen to bring the approximate solution into agreement 
with several numerically integrated flight path solutions. Substitution of equations (C5), 
(C6), and (C7) into equation (C4) yields the final equation for W p T( -,. 

The compressor and bypass ratios, P^/P^ and BPR, optimized for nuclear cruise 
were assumed constant during takeoff, climb, and acceleration, while a turbine temper- 
ature of 2160° R was selected. 


Descent Fuel W f ^ 

The descent fuel is based on a typical value of 2000 pounds per hour at descent con- 
ditions for a 300 pound per second airflow turbojet. Thus, for the more general case of 
a turbofan engine, 



2000 w a, T *d 
300 1 + BPR 3600 


(C9) 


The descent time t^ can be approximated by assuming that during descent the L/D 
ratio remains constant. The range covered during an idle power descent is then 
H d (L/D). If the velocity decreases linearly from Vq to the approach velocity V app , 

v = V 0-< V 0- V app>(^) < C10 > 


For large L/D, the range during descent is 


H d 




1 

2 


^0 + V app) 


(Cll) 


Thus, 




V 0 + V app 


(C12) 


Combining equations (C9) and (C12) yields 


• H, [-] 

W a T ^ l n 1 
W-ri j = 0. 0037 til yzl 

* ’ 1 + BPR V n + V, 


(C13) 


0 app 


Emergency Cruise Fuel Wp c 

To achieve a given range R at the design flight conditions, the required fuel is 


W F> c = D(sfc)t 

where the drag D is (W g - W F f ^,)/(L/D) and the time t is determined by 


(C14) 


t = 


* - «a(| 


(Cl 5) 


Hence, 


W 


F,c“ 



(sfc) 



(Cl 6) 


The sfc is calculated for the nuclear cruise cycle parameters P^/Pj, BPR, and T^. 


APPENDIX D 


AERODYNAMICS 

The variation of maximum fuselage width with gross weight was selected as 

B = 16. 4 + 11. 2(— i. - 0. 5) (Dl) 

\ 10 6 / 

The surface area of the fuselage group composed of the basic body, the reactor shield 
fairing, the gear sponsons, and the dorsal fin was given by 






1. 56 R 


sh 


+ 


B 

2 




W„ 

+ 750— i- 
10 6 


(D2) 


Examples of the variation with gross weight of fuselage width, length, surface area, 
and cargo bay length are given in table VII. The values for fuselage length to width ra- 
tios of 12 and 9 are shown as well as having the reactor assembly either within or exter- 
nal to the cargo bay. 


TABLE VH. - EXAMPLES OF FUSELAGE DIMENSION PARAMETERS 


w g’ 

lb 

B, 

ft 

Reactors exterior to cargo bay 
< R sh= 7 > 

Reactors interior to 
cargo bay (R gh = 0) 

S f> 

ft 2 

{If/ B= 12) 

l t> 

ft 

(Z f /B = 12) 

l c’ 

ft 

(Z { /B= 12) 

Sf, 

ft 2 

(Z f /B = 12) 

S f> 

ft 2 

(If/ B= 9) 

0. 5xl0 6 

16.4 

11 500 

196.7 

123 



1.0 

22.0 

21 100 

264.0 

165 

15 027 

13 640 

1.25 

24. 8 

26 750 

297. 5 

186 



1. 5 

27. 6 

32 500 

331.0 

207 
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The wing thickness ratio perpendicular to the quarter chord line was reduced from 
the value of 0. 18 for Mach numbers greater than 0. 65 according to the relation 


r w, 1,1/4 


0. 117 


Mq + - (M q - 0. 65) 
3 


(D3) 


The thickness ratio of the tail was maintained at 0. 12. 

The planform taper ratio of the wing was maintained at 0. 3 and the sweepback angle 
was calculated by equation (D4) which has been slightly modified from that given in ref- 
erence 11: 


A 


1/4" 


T>. 655 - 1. 77 t 


w,l, 1/4 t Q 


72 


cos 


-1 


1.794 


C L 


M, 



(D4) 


Table VIII gives examples of values of wing thickness ratio and sweepback angle with 
Mach number for a lift coefficient C L of 0. 6. 


TABLE VIII. - EXAMPLES OF WING PARAMETERS FOR LIFT 


COEFFICIENT OF 0. 6 


Flight 

Mach 

number, 

M q 

Wing thickness to chord 
ratio perpendicular to 
wing quarter chord line, 

T w, 1,1/4 

Wing sweepi'ack 
angle (quarter 
chord), 

A l/4> 

deg 

Sweepback correction 
to induced drag, 

0. 4 to 0. 65 

0. 18 

0 

1. 0 

0. 70 

. 151 

0 

1.0 

.75 

. 132 

22 

1. 06 

. 80 

. 117 

30 

1. 10 
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Figure 22 . - Parametric representation of engine nacelle and pylon strut. 


The surface area of the nacelle and pylon depicted parametrically in figure 22 was 


S = N ^ . f lE.. (bPR + 4. 5 + — ) 
w BPR + 1 \ n/ 


(D5) 


BPR = | -J“ 1 - 1 
D tip/ 


(D6) 


D 2 -A 

fan ' ttN 


F 

inst 


JL 49.4 A! A 


(D7) 


w. 


with F- f taken from equation (7). The friction drag of the nacelle and pylon was 

inst 

doubled to account for scrubbing and interference effects. 

The coefficient of drag due to lift was compiled as 
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where 


Q a = f(A) (See table VIII) 

Q x = f(X) = 0.0021 for X = 0. 3 
Q = 0. 006 

r 

The general equation used for the profile drag of major components was 


C D,o,i = A il- li I7 r - C r,o] C i 
l i 


(D10) 


The coefficient of friction Cj r was calculated using the method of reference 8 for com- 
pressible turbulent boundary layer as a function of Reynolds number based on the char- 
acteristic length of each major component. The coefficient C. is used to adjust the 
friction and roughness drag for other effects such as thickness ratio for the airfoil sur- 
faces (ref. 9): 


C 4 = 1 + 2t + 6 Or 4 


(Dll) 
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APPENDIX E 


STRUCTURAL WEIGHT ESTIMATION 

The gross weight at the beginning of cruise is used as the basis for structural weight 
calculations. This differs from the takeoff gross weight by the fuel consumed during 
takeoff, climb, and acceleration: 


W g = W g " W F, TC 


(El) 


The flight gross weight can also be expressed in terms used in calculating the weight of 
the structure as 


W g= W U +W w > a + W w < E2 > 

The term ^ conceptually includes everything attached to the fuselage other than the 
wing and can be further defined in terms of structural items, powerplant components, 
and payload as 


W f,Z= W f,a + W fuse + W t + W fix + W gear 


(E3) 


The total load available in the fuselage for the arrangement herein is the sum of the re- 
actor and shield .(including penetration penalty), W , and the payload, W T : 

cd J-J 


W f, a = W s+p+rea + W L 


(E4) 


The two terms remaining from equation (E2), W„, and W , represent the sum of the 
weight factors causing bending moment relief in the wing. Items other than the wing 
weight W,„ are included in W,„ „ which is defined, for the arrangement herein, as 

|V W y 

W w,a = W eng + w N +W F + tanks +W HXL < E5 > 


The procedures used for calculating the weights needed to evaluate equations (E2) to (E5) 
are now given or explained. The weights W s+fH _ rea , and W eng are known 

from separate routines such as appendixes A and B. The weight of the nacelles was es- 
timated using the engine-pod (nacelle-pylon) surface area from the configuration given in 
appendix D, and a weight of 3 pounds per square foot: 


W N = 3NS n 


(E6) 


The remaining item from equation (E5) that relieves wing bending moment is the chemi- 
cal fuel and its containment weight referred to herein as chemical fuel plus tanks: 

W F + tanks = >• 032 W F - W F, TC < E7 > 

The structural weight of the wing W w was calculated for X = 0. 3 using methods similar 
to those of reference 12: 


W = 
w 


1 S“„ n,b ( b c 

h. 86 S + 0. 00388 — + -i_® W’ 15. 6 + 3. 17 

i + n h / b\] w b 1q 6 g \ t 

Li (13 . 2 + 2.33-1] l s 10 L \ r ' 


? F+tanks (l2. 8 + 2. 06 - (W eng + Whxl + W„) L. 9 + 


where the structural span is 


cos A q 4 


the maximum wing root thickness is 


t = c T. 


r =c r T w, 1,1/4 cosA 1/4 


and the root chord is 


AR (1 + X) 


The structural or equipment weight items in equation (E3) were calculated as fol- 


lows: 


W fix = 10 000 + °- 03 W g 
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(E13) 


W = W 
gear "g 


/ W'^" 
0. 05 1 5 [ — I 
6 


V 


10" 


' Sui. S i. 

W t = 0. 012 [— + — jnjW 

S W S W; 


g 


(E14) 


The payload weight is given by 


W L = W cargo + W crew + W prov 


(El 5) 


The fuselage structural weight was calculated using a slightly modified relation from ref- 
erence 12: 


W fuse = c 3 I 3 +-V-Id? + C„n f d f (-i ||1.3+i 


d f / \ d f y 


f + 4 f u f 


2. 3 - k| + 2 X, /^W 


df/ '3 


17 4 j w M 


+ 5B 2 ^ - 4 A + 16. 9 B 2 + 10 [jt(B + df) + 8. 2 b] (E16) 


The last three terms in equation (E16) are simple-minded adjustments for cargo floor, 
pressure bulkheads, and front and rear loading door cutouts. 

The maximum depth or height of the fuselage was chosen as 


d f = 1.273 B + 1.25 R g 


Kf = K 


‘1/4 


c. g. 


(E17) 

(E18) 


K c. g. = °- 55 


C 3 - 3-2 



The airframe or structural weight estimated by this procedure corresponds to the usual 
empty weight minus the dry weight of the engines and is defined by 

W str = W w + w fuse + W t + W fix + W gear + W tanfcs < E19 > 
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APPENDIX F 


OPTIMIZATION PROCEDURE 

The optimization routine keeps track of the independent variables and the figure of 
merit. It perturbs the initially guessed independent variables individually, and stores a 
new value (or guess for the optimum value) for the independent variable if the figure of 
merit has been improved. After a pass through the entire array of independent variables, 
perturbation factors are changed and the process is repeated. When the change in the 
figure of merit at the end of the pass is less than some specified tolerance, the case is 
terminated and the results outputed (flow chart A, fig. 23). Of course, this procedure 
does not guarantee an absolute maximum. It may be local, depending on the choice of 
starting conditions, and several sets of starting conditions may be necessary to deter- 
mine the absolute maximum. Also, it is possible to obtain nearly the same optimum 
with different combinations of independent variables. This occurs since a narrow band 
of alternate selection exists within which the independent variables may be mutually 
traded without altering the figure of merit (within the convergence tolerance). 

The independent variables considered in this study are Pg/Pj, Pg/Pj, P4/P3, 

BPR, T^, T ex , T hx , W/A, Pjjx, P duc t> AR, and This choice is arbitrary as 

long as the independent variables plus fixed inputs are sufficient to completely specify 
the system. The choice of engine parameters is fairly typical and will determine engine 
performance if the shaft power extracted for pumping helium is known. In addition, the 
independent variables specify inlet and outlet air side conditions for the heat exchangers. 

Guessed values for the optimum independent variables including AR and with 
fixed choices of fuselage sizing parameters (appendix D) and an assumed value of 
yield an initial guess as to the drag of the aircraft. From this, the reactor power is 
calculated. Then, the reactor power together with T gx , T-^, W/A, and P Hg size the 
reactor and shield. These, in turn, adjust the fuselage and engine sizes and determine 
Cjy This Cq thus becomes a variable in a loop until the new assumed Cp agrees 
with the calculated one. 

Helium side temperatures and pressure drops plus air side conditions from the en- 
gine calculations are then used to size the heat exchangers and ducts. 

Finally, the numerical value of the power is needed to calculate the helium pressure 
drop through the reactor - hence, the shaft power to run the pump. Since the power is a 
function of the pressure drop and visa versa, another loop exists and must be iterated. 

The need for an overall optimization procedure becomes apparent when the engine 
air side, nuclear helium side, and airframe components are optimized independently. 

For example, if the engine parameters are optimized on the basis of thermal efficiency, 
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would increase beyond that obtainable from the heated helium. The limiting factor 
would have been the growth of heat exchanger weight and of the power required. If L/D 
is used to design the airframe, the structure weight is so high that the corresponding de- 
crease in engine size and reactor power still results in a net decrease in the value of the 
figure of merit. If helium side pressure drops are. minimized in order to reduce power, 
the heat exchanger and helium duct weights increase. These are just a few of the trade- 
offs involved. Engine parameters not only affect engine size, but also the power re- 
quired - consequently, the helium side conditions. These, in turn, affect reactor size 



(a) MAIN routine. 



(b) Flow chart A: optimization routine. 

Figure 23. - Flow diagrams for airplane-powerplant optimizing procedure. 
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and thus the fuselage dimensions, which affect the drag and hence the engine size, ad 
infinitum. 

Any other independent variable enters this same loop at some point and the need for 
a master program to oversee these effects becomes apparent. The overall process can 
be seen by examining the main, optimization, airplane design, and calculation of figure 
of merit flow charts (see fig. 23). 



(c) Flow chart B: design of airplane. 



(d) Flow chart C: calculation of figure of merit. 
Figure 23. - Concluded. 
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